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Abstract 

We investigate the phase diagram of antiferromagnetic spin ladders with ad- 
ditional exchange interactions on diagonal bonds by variational and numer- 
ical methods. These generalized spin ladders interpolate smoothly between 
the 5 = 1/2 chain with competing nn and nnn interactions, the S = 1/2 
chain with alternating exchange and the antiferromagnetic S = 1 chain. The 
Majumdar-Ghosh ground states are formulated as matrix product states and 
are shown to exhibit the same type of hidden order as the af S = 1 chain. 
Generalized matrix product states are used for a variational calculation of the 
ground state energy and the spin and string correlation functions. Numerical 
(Lanczos) calculations of the energies of the ground state and of the low-lying 
excited states are performed, and compare reasonably with the variational 
approach. Our results support the hypothesis that the dimer and Majumdar- 
Ghosh points are in the same phase as the af S = 1 chain. 
PACS numbers: 74.20Hi, 75.10Lp, 71.10+x 



1 



Typeset using REVTgX 



I. INTRODUCTION 



In the last few years spin systems consisting of a finite number n of interacting spin 
chains with S — | (now usually called spin ladders, consisting of legs and rungs, see fig. [I] 
for n = 2) have attracted considerable attention as recently reviewed by Dagotto and Rice!. 
^From the theoretical point of view, these systems are intermediate between one- and two- 
dimensional systems. Interest in these systems started with the two leg ladder! with isotropic 
antiferromagnetic couplings on the legs and rungs; in contrast to the 5* = \ — chain this 
ladder is characterized by a spin gap. The appearance of this gap is immediately clear in 
the limit of strong ferromagnetic interactions on the rungs (it then is identical to the gap 
characterizing the S = 1 antiferromagnetic (Haldane) chain) and in the limit of interactions 
on the rungs only (dimer limit); the gap, however, apparently persists not only for the 
isotropic antiferromagnetic ladderH but for nearly all values of the coupling constants 
except in the limiting case of two decoupled chains. In more recent investigations, af ladders 
with n legs have been shown to behave as spin liquids for n odd and as gapped systems for 
n evenliS. 

Experimental investigations have been performed on the n = 2 ladder substance 
(VO) 2 P20 7 i and also on the family of compounds Sr m _ 1 Cu m+1 2m ( m = 3,5,7...)! which 
approximately realize ladders with n — l(m + 1) legs. In these materials evidence for the 
existence of a spin gap has been found in susceptibility!'!, neutron scattering^ and magnetic 
resonance0lll experiments. The cuprates are of additional interest as candidates for high 
temperature superconductivity, which, however, remains undetected^. 

These experimental results have been successfully dealt with in a large number of theo- 
retical approaches using e.g. perturbation theory and exact numerical diagonalization00, 
mean field theory based on bond operators!! and the density matrix renormalization group 
approach!!!. 

Recently theoretical interest in spin ladder models with n = 2 has concentrated on the 
fact that these models allow to interpolate smoothly between seemingly very different spin 
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systems when the rung coupling is varied, such as between the S = 1 antiferromagnetic 
(Haldane) chain and the dimer state, i.e. models which are built of triplets, resp. singlets 
on the rungsQBEl. 

An interpolation between these two limiting models which avoids the singular point of 
two decoupled legs is possible when an additional interaction on one of the diagonal bonds is 
introduced as proposed by WhitJ and Chitra et aE3. This model is defined by the following 
Hamiltonian (see fig. [l]): 

H = H^ + ^H^ + (1+ 7 )H^ (1) 

= e s 1}j s 2sj hm = £ £ 4 s iJ+1 h& = £ s ld s w (2) 
j=i j=i 1=1 j=i 

Periodic boundary conditions are used and operators Sij denote spin-1/2 operators with the 
indices i — 1, 2 labeling the two legs and j — 1,2 ... L labeling the L rungs. When sites are 
labeled by a single index (running from 1 to 2L as also indicated in fig. [I]) it becomes clear 
that the spin ladder with additional diagonal couplings is a generalized model composed of 
the spin-1/2 chains with bond alternation and next nearest neighbour (nnn) interaction. 
This model includes in particular (compare with fig. 0): 

• the antiferromagnetic (af) S = \ Heisenberg chain with nnn exchange of relative 
strength ~ (a-axis, 7 = 0)), including in particular the Majumdar- Ghosh (MG) limitEH 
(a = 1,7 = 0) with two degenerate dimerized exact groundstates given as a product 
of singlets either on the rungs or on the diagonals, 

• the af S = I Heisenberg chain with alternating exchange (7-axis, a = 0), including in 
particular the dimerized chain (a = 0,7 = —1) with couplings on the rungs only and 
a product of singlets on the rungs as exact ground state, 

• the "regular" isotropic S = | ladder for 7 = — 1; the experimentally relevant cases 
mentioned above are expected to correspond to a = 2, i.e. to coupling of equal strength 
on the rungs and on the legs, 



• the Haldane chain with L sites for 7 — > —00, a > and nondiverging (leading to S = 1 
units with an effective af coupling |(1 + a) on the diagonals), 

• the isotropic S 1 = \ af Heisenberg chain with 2L sites (for a = 0, 7 = 0) and two 
decoupled af Heisenberg chains with L sites each (for a — > 00 with 7 finite). 

In addition the following lines in the phase diagram are of particular interest: 

• the line (a), 7 = 0, a > 0, is of particular interest since on this line a transition 
from the gap less isotropic af Heisenberg chain to the gapped MG chain occurs with 
increasing nnn coupling strength a. This transition has been investigated in detail@0 
(for a review see ref. |21| ) and the critical point has recently been located very accurately 
at \a cr ~ 0.2411670. Line (a) and its immediate neighbourhood (7 <C 1) is also of 
experimental relevance because of the discovery of the spin-Peierls compound CuGeO"3 
with alternating and strong nnn interactions^; 

• the straight line (b), which connects the dimer point to the MG point ((0,7) = 
(0, —1) — > (1, 0)); on this line one of the MG ground states, the state with singlets on 
the rungs remains an exact ground stateB 

In this paper we will discuss the properties of the generalized spin ladder in the space 
spanned by the interaction constants a > and 7. First we note that there is a symmetry 
transformation connecting the upper half-plane of the phase space of fig. ^| to the strip 
— 1 < 7 < 0: a translation of the upper leg by one lattice constant to the left leads to a ladder 
with the same symmetry, but different coupling constants. The symmetry transformation 
exchanges the roles of rungs and diagonals; apart from a change in energy scale by a factor 
1+7 (which is qualitatively irrelevant for 7 > —1), it is expressed in the following way in 
terms of the parameters of the Hamiltonian: 

1+7 1+7 
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The symmetry line of this transformation is the a— axis 7 = 0, points on this axis are 
mapped onto themselves. Owing to the presence of this symmetry we will consider in the 
following the region 7 < of the a — 7— plane only. The strip — 1 < 7 < is related to 
the upper half-plane by the symmetry of eq.(^), the section 7 > — 1 of the a — 7-plane 
has af coupling and is identical to the phase space discussed by Shastry and Sutherland^. 
Including the region 7 < — 1 adds the possibility of ferromagnetic coupling between the two 
legs, this possibility has so far been discussed for special cases onlyi!0. 

In the following we will present results for the generalized spin ladder as defined by the 
Hamiltonian of eq.(|l]) from a variational approach based on the concept of hidden order. 
These results will be supplemented by and compared to the results from exact diagonaliza- 
tions for finite ladders with up to 14 rungs. The question of hidden order in generalized 
spin ladders comes up naturally since the af S = 1 chain, the prototype system where 
hidden order has emerged first, is realized in the limit 7 — > — 00. In the Haldane chain@ 
hidden order has been made transparent by Kennedy and Tasaki0, using a nonlocal uni- 
tary transformation. Ground states with complete hidden order were then written down 
in the form of matrix product (MP) ground statesil and used as basis for variational cal- 
culations. The new feature of ladders is that singlets on the rungs get mixed in when one 
moves from the Haldane limit 7 — > —00 towards finite values of 7. A generalization of the 
Kennedy- Tasaki-transformation0 was proposed by Takada and WatanabeH for the isotropic 
spin ladder (7 = —1, varying a). Explicit definitions of hidden order in spin ladders were 
recently discussed by Nishiyama et al0 and White!. 

Our variational approach starts from the observation that the two exact ground states 
of the Majumdar-Ghosh Hamiltonian can be written in the form of matrix product states. 
We then generalize the matrix product state to allow the ansatz to include the approximate 
ground state of the Haldane chain as given by Affleck et a and written as matrix product 
state by Kliimper et afi Continuous variation of the parameters then smoothly connects 
the seemingly unrelated limiting cases of the Majumdar-Ghosh and dimer points on the one 
hand and of the Haldane chain on the other hand. We expect this ansatz to make sense as 



a frame for a qualitative description in a large part of the half plane a > 0. 

In section 2 we will introduce the variational matrix product states defining the hidden 
order in the analogous way as for an S=l chain but replacing the contribution of the S z = 
component of the rung triplet state by a linear combination of this component with the 
rung singlet. The exactly known ground states of the Majumdar-Ghosh chain are particular 
examples of this ansatz. This definition of matrix product states imposes a particular type 
of hidden order which we will describe. In section 3 we present the results of variational 
calculations for the ground state energy and for ground state correlation functions obtained 
by using this scheme. In addition the string correlation function as proposed in ref. [5] will 
be computed. In section 4 we present numerical results for the ground state energy and 
the gap along various paths in the phase diagram (shown as dashed lines in fig. |2|). These 
results are obtained from exact diagonalization of finite ladders with up to 14 rungs and 
periodic boundary conditions using the Lanczos algorithm; they will be discussed in compar- 
ison to the results of the variational calculations. Our results imply that a smooth variation 
of parameters leads one from the Haldane-phase to the dimer-phase; they therefore sup- 
port the hypothesis that these two limiting models are not separated by a phase transition 
and that the only critical behaviour in the phase diagram considered is on the gapless line 
7 = 0, < a < a cr . Special attention will be paid to the neighborhood of this line in the 
numerical calculations. Conclusions will be given in section 5. 



II. MATRIX-PRODUCT STATES FOR SPIN LADDERS 

We start by a discussion of the Majumdar-Ghosh chain, i.e. a = 1,7 = 0. For periodic 
boundary conditions the two ground states are known exactly: they are dimerized config- 
urations with singlets on either the rungs, |0 r ), or the diagonal bonds, 1 0^) . If we use the 
following representation of states on the j-th rung, 

l*)i = ^(IU>i-UT>i), 
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\toh = ^(| Uh + I IT),), lUh = I TT)j> = I U>;, (4) 

the two ground states can be writen in form of MP-states 

|0r> = ni*>i \O d ) = (-l) L Tr(flgA (5) 
j'=i \i=i / 

with the matrices ^ defined by 

1 ( \s)j + \to)j -V2\t+)j \ 
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(6) 



' \ V5\t-)j \s)j - \k)j j 
In the representation of eq.(^jj) both MG ground states are given in terms of singlets resp. 
triplets on the rungs; this appears complicated for the state \0d) but actually it serves to 
realize the presence of hidden order in the MG ground states: The matrix representation 
of the ground state according to eqs.([||)]) has components of the following string structure 
only 

...t+ t- t+ t- (t +s)(t +s) (t +s)(t +s) t+ t- t+ t- t+ (to-s)(t -s) (to-s)(*o-») t- t+ t- t+..., (7) 

i.e. the sequence . . . i+ t_ t+ t- ... is interrupted either by an arbitrary number of sites with 
| to + s) following a site with or by an arbitrary number of sites with |i — s) following 
a site with \t + ). 

The symmetry transformation of eq. (|^) transforms the ladder into itself on the line 7 = 
and the two states |0 S ) and 1 0^) remain degenerate for a / 1; they are, however, ground 
states only at the Majumdar- Ghosh point. When moving off the axis 7 = this symmetry 
is destroyed and at most one of these states survives as ground state. In particular, on 
the line (b), (0,7) = (0,-1) — > (1,0) the ground state is |0 r ), i.e. composed of rung 
singlets. This 'disorder line' is therefore characterized by a vanishing of the correlation 
length. Under the symmetry transformation (|3|) this line (b) transforms into the line (b') 
(a, 7) = (1,0) — > (l,oo) with a dimerized ground state with singlets located on diagonal 
bonds. As discussed in the introduction it is sufficient to discuss only the region 7 < in 
the phase diagram of fig. |2|. 



When we move away from the MG point, the state |0 r ) has no freedom to adapt but 
the state \0d) can be used as the basis for a variational calculation by allowing different 
amplitudes for the singlet and the triplet contributions. Rotational invariance requires the 
following form for the matrix gj 
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Tr ( n 93 ] > m 

u=i 



f a\t )j + b\s)j —aV2\t+)j ^ 



(8) 



v aV2|t-)j -(a|t )i - b\s)j) J 
This follows from the fact that the three states 

IO = -^(1 tt> - I U», IO = ^(ltt> + IU», 

r) = ^(lTI) + IIT))=:|M (9) 

form a vector and can only appear in the 2x2 matrix gj in the combination J2 a & a \t a ) in the 
case of rotational invariance (a a are the Pauli spin matrices). In eq.(f|), a and b are complex 
variational parameters, subject to the normalization condition 3|a| 2 + \b\ 2 = 1. Thus the 
parameters entering the variational calculation are one amplitude and one phase. 

The ansatz of eq.@ can be shown to reduce to the state introduced by Takada and 
WatanabelH when their state is written in original spin space (no unitary transformation) 
and perfect generalized hidden order is introduced. If we put b = in eq. (§) we obtain the 
AKLT stateB the exact ground state of the af S — 1 chain with biquadratic exchange of 
strength ~. 

It has to be noticed that the variational ansatz of eq.(G|) may be used with the matrices 
gj defined either on the rungs or on the diagonals and as a first step of the variational 
calculation it is necessary to find out which of these possibilities leads to lower energy. It 
turns out that for 7 < (which is the region we will consider) lower energy is obtained when 
gj is defined on the diagonals, whereas defining gj on rungs gives the lower energy for 7 > 0. 
On the symmetry line 7 = eq.(§) is an approximate description of the two equivalent 
groundstates for a cr < o -C 00. The ansatz of eq.(|) allows to reproduce the exact ground 
state on line (b) (7 = — 1 + a), the variational approach is therefore exact on this line. 



On the one hand, the wave function of eq.([D is rather general since it can smoothly 
interpolate between dimerized states and AKLT-states, i.e. states which correspond to points 
which are far separated in the phase diagram. On the other hand one cannot expect this 
approach to work for the whole a-7-plane since it is manifestly inadequate for special points 
like the ones corresponding to the HAF or to two decoupled chains, which have gapless 
spectra and power-law decay of correlations functions. 

We finally note that there is another possibility of constructing a rotationally invariant 
wave function of the MP type which is related to the resonating valence bond structure also 
discussed for spin ladders: a MP wavefunction with the following alternating structure 
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RVB 



) = Tr 



II 9tk~i92k > af 



1 



k=l 



\ 



(10) 



±\s), + \t ) j -V2\t+)j 
v^lt-), ±\s)j - |t ). 

i.e. the previous structure with diagonal elements interchanged in every second g-matrix, 
describes a state with singlets located on the legs. Thus our ansatz is also connected to RVB 
states on ladders as introduced in refs. B2,B3. 



III. VARIATIONAL CALCULATIONS FOR GENERALIZED SPIN LADDERS 

In this section we present the calculation of the energy, of the spin correlation function 
and of the string (hidden order) correlation function using the variational state of eq.(§). 
We discuss results for these quantities in the a — 7-plane and compare to the results of alter- 
native approximate approaches for special points. A comparison to the results of numerical 
calculations will be given in section 4. 

The calculation of the ground state energy and of the ground state correlation functions 
can be done in complete analogy to the corresponding calculations for the S = 1 af chainil. 
The variational energy, i.e. the expectation value of the Hamiltonian of eq.(|l]) with the 
wavefunction of eq.(H) and the matrix gj defined on diagonal bonds, is obtained as 

E^L = _ 3{y4 4 + 2A 3 B ^ + A 2 B 2 cog 2 ^ 
Lj 
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-3a{A 4 - A 2 B 2 cos 2 <p} + ~(1 + 7 ){A 2 - B 2 }, (11) 
A = | a |, B = \b\, 

with the norm 3 A 2 + 5 2 = 1 as an additional constraint, if is the relative phase of the 
original complex amplitudes a and b. 

E var reduces to the expression given by WatanabeEl for 7 = —1 and to the expression 
given by Takada and WatanabeH for 7 — > —00, a oc (1 + 7) for perfect generalized hidden 
order. The procedure of minimizing E var is performed easiest by first minimizing with 
respect to cosy? and then with respect to the ratio u = B/A. The analysis shows that the 
absolute minimum always corresponds to cos ip = ±1, and the amplitudes a and b can always 
be taken real; the remaining equation for u is of third order. Thus the minimum can be 
calculated exactly for general values of a and 7. 

The simple variational wave function of eq.(||) also allows the calculation of the ground 
state correlation functions. They exhibit the exponential decay characteristic for MP states. 
Owing to isotropy there is only one correlation length which is obtained as 

In order to discuss the hidden order for the ladder system, we use the definition given in 
ref. U, which reproduces the known results in the S = 1— chain limit. Using this definition, 
we obtain the following result for the string correlation function (in the thermodynamic limit 
L — > oo) 

1*7(01 := KiSlo + SlAe^>l k+ s lk+l)(Sh + = 4 (J _ ^ 2)2 _ (13) 

In fig. [3] we show results obtained from the MP state ansatz for points on the lines 
7 = 0, —1. In the MG-case (7 = 0, a = 1) and at the dimer point the variational ground 
state is exact, we have A = B = |, cosy? = 1. The wavefunction is identical to that of 
eqs.(|]|l), one of the two exact MG ground state wavefunctions with energy E MG = —\L. 
The correlation lengths £mg — idimer vanish. Figs. ||(a,b) show the ground state energy 
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and the correlation length on the lines 7 = 0,-1. Discussing our results on the line 7 = 
i.e. for the general af chain with nnn interactions, we have to restrict the discussion to the 
range a > a cr ~ 0.5 since the MP ansatz is no more appropriate when the gapless system 
is approached. Nevertheless we note that the variational ground state energy for a = a cr 
is Eq/L = —0.79354, i.e. within 1% of the exact value0. In the region to the right of the 
disorder line the relative error increases to typically 5%. The correlation lengths remain 
very small as is typical for MP wave functions. 

The string order parameter is determined by the singlet weight B 2 which is shown in 
fig.^c, we have B 2 > | to the left and B 2 < \ to the right of the disorder line (b). The string 
order parameter for the two degenerate MG states takes the values go r {l) = 1/4, 9o d (l) = 
(obtained from B Qr = 1/2, B 0d = 1). The fact that the string correlation is different for the 
two equivalent ground states is related to the asymmetric definition of the string correlation 



with respect to rungs and diagonals in eq.(|l"3"|). 



For the symmetry line 7 = 0, another variational wavefunction of RVB type (which is 
also exact at the MG point) was proposed by Zen and Parkinsonil. When the results for 
the ground state energy are compared for these two approaches, for a > 1 lower energies 
are found from the RVB approach and for a < 1 from the present MP state approach. The 
correlation length for a > 1 is similar in the RVB and the MP state approaches, wheras 
for a < 1 the MP ansatz results in a larger increase in £. The behaviour of the correlation 
lengths in this region of the phase diagram might be influenced by the possible emergence 
of a spatially modulated ground state (spiral phase) at a = 100. The MP state ansatz 
cannot by construction reproduce such a behaviour and the alternating ansatz as described 
above in eq. (|T0|) does not lead to lower variational energies either (except for large values of 
a). A more detailed discussion of this aspect will be given when we discuss our numerical 
results in section 4. 

We now turn to a discussion of results for the regular ladder system (7 = — 1, a = 2). 
From the variational approach we obtain Eq/L « —1.102 and £ ~ 0.815 for this isotropic 
ladder. These results may be compared to those for a variational RVB wave function which 
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was proposed by Fan and Maa and leads to a ground state energy of Eq/L ~ —1.112 and a 
correlation length of £ ~ 0.238. The difference between the two variational energies of about 
1% appears negligible when it is compared to the exact (numerical) result Eq/L « —1.156. 
Both variational methods lead to correlation lengths considerably smaller than £ « 3.19 as 
computed in ref. This is consistent with the well-known fact that in approaches using 
MP states the correlation length is notoriously underestimated (compare the AKLT value 
of £ ~ 0.9102 to the exact correlation length £ « 6.2 in the af S = 1 chain). 

Finally we want to illustrate in the set of figs. §](a-c) how the limit of an antiferromagnetic 
5 = 1 chain is included in the present scheme for 7 — > — 00. In this limit it is very unfavorable 
to have singlets on the diagonal bonds, so minimization leads to B — > 0. We show the ground 
state energy and the singlet weight with increasing | — ^ | for a = 2 in fig. |](a) an d in fig. £|(c) 
respectively. Combining these graphs with the corresponding ones from fig. |3| we can follow 
a continuous path from the dimer point to the Haldane limit. The smooth variation of the 
physical quantities on such a path illustrates the similarity of the Haldane phase and the 
dimer phase in the present approach. 

The minimum of the variational energy in the limit appropriate for the af S = 1— chain 
has the following form 

^%!-H^-^+^y l (ii+ ^ i>>i) 

82 » 3(TW (14) 
The first term in E$ =1 is just the internal energy of the triplets representing the internal 
energy of spins S — 1 at the L sites. We see that minimization in the limit of infinite 
negative 7 leads to the AKLT-state with 

A 2 = \, E = -p eff , J e// = (a + l)/4. (15) 

This is the identical result as obtained from direct variational calculations for the S — 1 
chain0. 
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For the S = 1— chain the hidden order is characterized by the string order parameter as 



introduced in Ref. |35|. In our approach, the AKLT value g(oo) — | is reproduced in the limit 



of zero singlet weight, i.e. for B = 0. As is seen from the monotonic behaviour of the singlet 



weight in fig. 0(c) in combination with eq.(O), the string order parameter g(oo) decreases 



monotonically when 7 is increased from —00 towards 0. 

IV. NUMERICAL RESULTS 

We have computed the ground state energy per rung and the energy gap for different 
values of the nnn interaction a and the alternating bond exchange 7 in order to compare 
these exact results to those obtained in the variational calculation of the last section in the 
various regions of the phase diagram. The paths which we have investigated numerically are 
shown in fig. ||] as dashed lines. The results for points in the upper half plane follow from 
the symmetry property given in eq.(||). 

Numerical calculations were done on the MPP CRAY T3D SC256 of the Zuse Computing 
Center Berlin. We used the Lanczos technique to determine the ground state energy per 
rung as well as the singlet-triplet energy gap on 2 x L lattices for L — 6,8, 10, 12, 14. In 



order to extrapolate the results for finite L to the bulk limit we used the ansatz of ref. [bl 



f{L) = f (00) + c e- L / L °L-v, (16) 



i.e. a power law in L multiplied by an exponential. Following ref. [T4| the Lanczos data were 
fitted using p = 2 for the ground state energy per rung and p = 1 for the gap energy. The 
typical relative mean square error of these fits was about 10~ 10 for the ground state energy 
and about 10~ 7 for the gap energy. The fit parameter Lq reflects the characteristic length 
of the system and corresponds to the correlation length in the chain under consideration. 
We have checked the accuracy of our program by recalculating data for the isotropic ladder 



(a = 2, 7 = —1) and have fully reproduced the results of ref. [14. 

In our numerical calculations we have concentrated on two different aspects of generalized 
spin ladders: 

13 



• On the variation of the ground state energy and of the excitation gap when we move 
from the dimer state to the Haldane limit. These data allow to estimate quantitatively 
the accuracy of the variational approach based on MP states and to discuss the issue 
whether these two limiting states are separated by a phase boundary. 

• On the neighbourhood of the gapless line (a) connecting the two conformal points 
7 = 0, a = (HAFM) and 7 = 0, a = a cri ~ 0.5. Such data are expected to 
contribute to an understanding of the crossover from the gapless phase on the critical 
line 7 = 0, a < a cr to the gapped phase(s) of the MG/dimer- and Haldane-type. 

We present the variations of the ground state energy, the gap energy and the length Lq 
from the dimer point to the Haldane limit on the paths shown as dashed lines in fig. |2|. 
Our numerical results for various chain lengths and their extrapolation to the bulk limit 
are given in tables | and |J and displayed in figs. |3] and |j. An extrapolation to the limit 
7 — > —00 was performed with a fourth order polynomial fit to E /L — (1 + 7)/4 for the 
points a = 2, 7 = —2.2, —2.6, —3.2, —4.4, —13. This extrapolation resulted in 



with J e ff = (a + l)/4 = 3/4. Both numbers are in excellent agreement with the well-known 
results for the ground state energy and the gap of the Haldane chainS. The analogous 
extrapolation for the characteristic length L (see eq.(|l6])) is shown in fig. f|(d); L increases 
monotonically when we approach the Haldane limit on our path. The variation of L Q corre- 
sponds to the variation of the correlation length £ although the limiting value of L is found 
somewhat larger than the accepted numerical value for the correlation length of the Haldane 
chain £ « 6.2. 

The numerical results are in agreement with the observation already made for the varia- 
tional results: On the continuous path through the phase diagram the quantities E , A, L 
vary smoothly and there is no indication of a phase boundary in the variation of these quan- 
tities. These data therefore support the hypothesis that the dimer point, the Majumdar- 
Ghosh point and the Haldane limit all are in the same phase. The gap characterizing these 




1 



A 



0.4106 



(17) 
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three limiting chains as well as the intermediate systems, including the regular spin ladder 
thus appears only quantitatively different and these systems therefore should be considered 
as manifestation of basically the same quantum condensation phenomenon. 

We now turn to the second region of interest in the phase diagram, the area around 
the line of vanishing gap connecting the two known conformal points. We calculated the 
energies of the ground state and of the first excited states 

• for 0.3 < a < 0.6 and small values of 7, allowing an extrapolation to 7 = 

• for values on the line 7 = —a, i.e. on a path from the gapless HAF towards the gapped 
phase at larger values of a, 7. 



In fig. |5] and in table |T| we show the behaviour of the gap when the gapless line, 7 = 0, is 
approached, i.e. the vanishing of the gap on the critical line and the (small) finite gap energy 
for a > a cr . For |-y| 1 and two values of a < a cr we have analyzed the results assuming 
a power law behaviour, A oc \5\ y , (where 5 = 7/(2 + 7) is the dimerization parameter) in 
order to compare with the prediction y = 2/3 by Cross and FisherS From a log-log plot 
(see insert in fig. |5p we find y ~ 0.862 for a = 0.30 and y ~ 0.795 for a = 0.44, i.e. an 



a-dependent exponent which is somewhat larger than predicted in ref. |37. 

The gap close to the conformal line is found at wavevector k = 0, whereas the gap in the 
Haldane phase is known to occur at wavevector k = tt (the wavevector k is defined through 
the factor exp(tk) which multiplies the wavefunction upon the transformation j — > j+1, i.e. a 
shift of rungs by one unit). To discuss this crossover in the wavevector of the lowest excitation 
we have calculated the lowest excitation energies for these two wavevectors, Ak=o, Afc =7r 
along the path 7 = —a. The results of an extrapolation to the thermodynamic limit as 
described above are shown in fig. || It is seen that the quantity Afc =0 — A^ =7r changes sign 
for a ~ 0.55 i.e. close to the line 7 = — 1 + a which connects the dimer and MG points. 
These results for the gap (and the corresponding results for the ground state energy) join 
smoothly to the results presented before when values a = 2 are reached. Actually the 
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excitation energies on the disorderline 7 = — 1 + a can be calculated in perturbation theory 
in a to a high degree of accuracy for a C 1. For k = it we have the exact result Afc =7r = 1 
on the whole disorder line, as will be published elsewhere. 

Generally, our numerical results confirm the variational results in that variations over 
most of the phase diagram are smooth and do not indicate any phase boundaries. However, 
we find that the excitation spectrum does change qualitatively along a path from the critical 
line to the gapped phase; this may be another aspect of the fact that the maximum of the 
structure factor S(k) is found at finite k and may be related to the existence of a spiral phase 
as speculated beforeil'0. In order to clarify this point, more detailed results on the excitation 
spectrum in this region of the phase diagram should be obtained. However, since a reliable 
finite size analysis is not possible for general values of the wavevector we cannot discuss 
the character of the complete excitation spectrum with comparable accuracy and therefore 
cannot decide whether the minimum excitation energy occurs at some finite value of the 
wavevector. In order to provide some first information we present in fig. [7] the excitation 
spectrum for the finite ladder with 2 x 12 sites at the points a = —7 = 0.54, 0.56, 0.58. 
The minimum excitation energy is found at k = for a = —7 = 0.54, at k = tt for 
a = —7 = 0.56. For a = —7 = 0.56 the minimum is found at finite wavevector, k « 27r/3, 
although dispersion is very weak and it is not clear whether the minimum at finite wavevector 
will survive an appropriate extrapolation to the infinite system. Actually, this result may be 
more accurate than one might expect from the small size of the system since we found in the 
finite size extrapolations for k = 0, tt that the results for the finite 2x12 site system differ 
very little (relative difference 10 -6 ) from the thermodynamic limit. For a full understanding 
more calculations are clearly required. 



V. CONCLUSIONS 

We have investigated the phase diagram of a generalized spin ladder with additional 
interaction on diagonal bonds. This model includes several known cases, in particular the 
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S = 1 chain in the Haldane phase and the dimer and Majumdar-Ghosh chains belonging to 
a dimerized phase. This generalized spin ladder was studied with a variational wave ansatz 
and numerical techniques. Using the concept of hidden order throughout the whole phase 
diagram (with the exception of the line of gapless points), a rotationally invariant matrix- 
product state was proposed for the ground state of these ladders. The ground state energy, 
and spin and string correlations were calculated. The comparison with numerical results 
for the ground state energies leads to agreement typically within about 4 - 5%, consistent 
with the qualitative nature of our approximation. In the S — 1 chain limit the well known 
AKLT state and the variational energy of 4/3 per site are reproduced and extrapolation 
of the numerical results leads to excellent agreement with the known values for the ground 
state energy and the gap of the Haldane chain. A string correlation function was chosen 
such as to reproduce in the corresponding limit the AKLT value of 4/9 for the factorizing 
approximation to the S = 1 chain. This string correlation function is directly related to the 
singlet weight in the MP wave function and decreases monotonically when one moves from 
the Haldane limit towards dimerized states. 

^From both the numerical and the variational calculations we find on paths connecting 
the dimer and Majumdar-Ghosh points to the Haldane limit smooth variations of all quan- 
tities considered. We are therefore lead to conclude that these dimerized af S = | chains 
and the af S = 1 chain are in the same phase and that the gap in the excitation spectrum of 
these chains is of the same nature. On the other hand, our numerical data for the excitation 
spectrum close to the line connecting the dimer and the MG points show a shift of the 
wavevector of the minimum excitation energy from k = to k = it as well as an indication 
that the minimum may be found at intermediate wavevectors. Thus the question of the 
existence of a spiral phase remains unresolved and requires further investigation. 

Using the present approach elementary excitations can be constructed as soliton-like 
states as for dimer-like systems^ and for the S = 1 chainS. Our approach can also be 
extended to include anisotropic couplings; in particular, similar as in the af S — 1 chain, 
there should be a critical value for the Ising-like anisotropy on the legs, at which the system 
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exhibits long range order even for 7 — > —00. These further applications of the MP approach 
are now under investigation. 
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TABLES 

TABLE I. Ground state energy Eq and gap energy Afc =7r for 7 = —1 and various values of a. 
TABLE II. Ground state energy Eq and gap energy Afc =7r for a = 2 and various values of 7. 
TABLE III. Gap energy A^ = o close to the gapless line (a). 
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FIGURES 



FIG. 1. Structure of the generalized spin ladder with additional diagonal interaction. 

FIG. 2. The phase diagram of the generalized spin ladder. Paths where numerical calculations 
have been performed are shown as dashed lines. 

FIG. 3. (a) Ground state energies, (b) correlation lengths and (c) singlet weight on the lines 
7 = and 7 = —1. Variational approach (full lines) vs. numerical results (O). 

FIG. 4. Approach to the Haldane limit on the line a = 2: (a) ground state energy, (b) gap 
energy, (c) singlet weight, (d) characteristic length Lq. 

FIG. 5. 7-dependence of the energy gap Afc = o close to the gapless line (a): 
a = 0.30(+), 0.44(x), 0.52(0), 0.60(D). Insert: Scaling of the gap with the dimerization pa- 
rameter 5 = 7/(2 + 7) for a = 0.30, 0.44. 

FIG. 6. Energy gaps Afc = o and Ak= w along the line 7 = —a. 5 = (a 2 + 7 2 ) 1 / 2 . Full lines are 
guides to the eye. 

FIG. 7. Excitation spectra e(k) close to the disorder line, —7 = a = 0.54, 0.56, 0.58 for L = 12. 
Insert: Expanded scale for —7 = a = 0.56. Full lines are guides to the eye. 
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Table 1 



Ground state energy 


7 = 


AT = 12 


N = 16 


N = 20 


N = 24 


N = 28 


N oo 


-1.0 


-1.168874 


-1.160406 


-1.157719 


-1.156744 


-1.156363 


-0.578020 


-2.2 


-1.442089 


-1.431225 


-1.427342 


-1.425759 


-1.425061 


-0.712146 


-2.6 


-1.536390 


-1.525068 


-1.520950 


-1.519238 


-1.518469 


-0.758797 


-3.2 


-1.679514 


-1.667673 


-1.663286 


-1.661423 


-1.660568 


-0.829772 


-4.4 


-1.969666 


-1.957152 


-1.952409 


-1.950343 


-1.949366 


-0.974069 


-13.0 


-4.095278 


-4.081498 


-4.076053 


-4.073567 


-4.072326 


-2.035283 


7 — y — oo 




-(7 + l)/4 - 0.52759 


Gap 


-1.0 


0.626569 


0.557398 


0.528107 


0.514784 


0.508496 


0.501711 


-2.2 


0.605725 


0.527135 


0.489590 


0.470340 


0.460095 


0.445421 


-2.6 


0.600401 


0.519686 


0.480410 


0.459877 


0.448724 


0.431682 


-3.2 


0.593758 


0.510573 


0.469294 


0.447250 


0.435003 


0.415168 


-4.4 


0.583965 


0.497512 


0.453582 


0.429489 


0.415715 


0.391956 


-13.0 


0.559096 


0.466436 


0.417295 


0.389032 


0.371958 


0.337796 


7 — > — oo 




0.30797 



Table 2 



Ground state energy 


a = 


N = 12 


N = 16 


N = 20 


N = 24 


N = 28 


N oo 


0.0 




0.75 


0.4 


0.766403 


0.766398 


0.766398 


0.7663978 


0.7663978 


0.7663978 


0.8 


0.819531 


0.819281 


0.819256 


0.819254 


0.819253 


0.819253 


1.2 


0.909511 


0.907827 


0.907510 


0.907443 


0.907428 


0.907424 


1.6 


1.029126 


1.024466 


1.023234 


1.022867 


1.022749 


1.022687 


Gap 


0.0 




1.0 


0.4 


0.821753 


0.821709 


0.821707 


0.821707 


0.821707 


0.821707 


0.8 


0.692602 


0.690373 


0.690074 


0.690032 


0.690026 


0.690025 


1.2 


0.619862 


0.605463 


0.601817 


0.600848 


0.600582 


0.600482 


1.6 


0.602390 


0.563670 


0.549962 


0.544825 


0.542838 


0.541330 



Table 3 



Gap energy 


a = 


7 = 0.10 


7 = 0.05 


7 = 0.03 


7 = 0.015 


7 = 0.00 


0.30 


0.2511807 


0.1497976 


0.0971629 


0.0526672 


< 0.00003 


0.44 


0.2799086 


0.1720255 


0.1167782 


0.0660249 


< 0.00003 


0.52 


0.3002298 


0.1864918 


0.1310506 




0.0009896 


0.6 


0.3251090 


0.2082006 


0.1507519 




0.0023626 
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